Optical forces arising from phase gradients 
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We demonstrate both theoretically and experimentally that gradients in the phase of a light field 
exert forces on illuminated objects, including forces transverse to the direction of propagation. This 
effect generalizes the notion of the photon orbital angular momentum carried by helical beams of 
light. We further demonstrate that these forces generally violate conservation of energy, and briefly 
discuss some ramifications of their non-conservativity. 



Light's ability to exert forces has been recognized since 
Kepler's De Cometis of 1619 described the deflection of 
comet tails by the sun's rays. Maxwell demonstrated that 
the momentum flux in a beam of light is proportional to 
the intensity and can be transferred to illuminated ob- 
jects, resulting in radiation pressure that pushes objects 
along the direction of propagation. This axial force has 
been distinguished from the transverse torque exerted by 
helical beams of light carrying orbital angular momen- 
tum (OAM) 0. We demonstrate theoretically and con- 
firm experimentally that OAM-induced torque is a spe- 
cial case of a general class of forces arising from phase 
gradients in beams of light. We also demonstrate that 
phase-gradient forces are generically non- conservative, 
and combine them with the conservative forces exerted 
by intensity gradients to create novel optical traps with 
structured force profiles. 

Our experimental demonstrations of phase-gradient 
forces make use of extended optical traps created through 
shape-phase holography [2U Sljj in an optimized [5] holo- 
graphic optical trapping [6|, system. Holographically 
sculpted intensity gradients enable these generalized op- 
tical tweezers |8[ to confine micrometer-scale colloidal 
particles to one-dimensional curves embedded in three 
dimensions. Independent control over the intensity and 
phase profiles along the curve then provide an ideal model 
system for characterizing the forces generated by phase 
gradients in beams of light. 



OPTICAL FORCES DUE TO PHASE 
GRADIENTS 

The vector potential describing a beam of light of fre- 
quency uj and polarization e may be written in the form 



A(r, i) = A(r) exp (i uot) e. 



(i) 



Assuming uniform polarization (and therefore a form of 
the paraxial approximation), the spatial dependence, 



A(r) = u{r) exp (i$(r)) . 



(2) 



is characterized by a non- negative real- valued amplitude, 
u(r), and a real- valued phase 3>(r). For a plane wave 
propagating in the z direction, <£(r) = — kz, where 



k = n m u;/c is the light's wavenumber, c is the speed 
of light in vacuum, and n m is the refractive index of the 
medium. Imposing a transverse phase profile <p(r) on the 
wavefronts of such a beam yields the more general form 



$(r) = -k z (r)z + <p(r), 



(3) 



where z • V<£ 
wave vector, 



and where the the direction of the 



k(r) = k z (r) z + V<p, 



(4) 



now varies with position. The associated electric and 
magnetic fields are given in the Lorenz gauge by 



E(r,t) 



dA(r,t) 
dt 



and 



H(r,t) = - V x A{r,t), 



(5) 
(6) 



where [i is the magnetic permeability of the medium, 
which we assume to be linear and isotropic. Following 
the commonly accepted Abraham formulation [9], the 
momentum flux carried by the beam is 



g (r) = -$1{E* x H} 



k 



fin n 



J(r)V$, 



(7) 

(8) 



where I(r) = u 2 (r) is the light's intensity, and where 
we have employed the gauge condition V • A = 0. This 
result is somewhat controversial [lOj despite a century of 
study, with the alternative Minkowski formulation differ- 
ing from Eq. (|8]) by a factor of n^. Both formulations 
are believed to yield consistent results for the force on an 
illuminated object, but differ in their predictions for the 
recoil of the medium [9[ . 

Applying Eq. (j3]), g(r) separates naturally into an axial 
component g z (r) = [kl(r) k z j (/m m )] z and an additional 
component 



9±(r) 



Hn Tl 



■ J(r) V<p 



(9) 



transverse to the optical axis. For small objects, the 
transverse force arising from g_\_(r) may be approximated 

as 



f ± (r) =<r(a)g ± (r), 



(10) 



2 



where a (a) is the cross-section describing the light-matter 
interaction. For a dielectric sphere of radius a, [H, [ll|, 



cr(a) = — /c 4 a b . 

3 V m ■ 



(ii) 



where m = n p /n m is the ratio of the particle's refractive 
index, n p , to the medium's, n m . 

More than a decade ago, Allen, et al [l| pointed out 
that the helical phase profile, (pe(r) = 10, imbues a beam 
of light with an orbital angular momentum flux, r x g ± , 
amounting to £h per photon. Here, 6 is the azimuthal 
angle around the optical axis, and £ is an integer describ- 
ing the wavefronts' helical pitch. This orbital angular 
momentum is distinct from the photons' intrinsic spin 
angular momentum p], EH, HI, 14, Ti|, 16j]. Through it, 
even linearly polarized light can exert a torque around 
the optical axis [13, EH, El[ • Equations ([9]) and ([TO]) re- 
veal the torque exerted by helical wavefronts to be a par- 
ticular manifestation of the more general class of trans- 
verse forces arising when phase gradients redirect radia- 
tion pressure. 

Intensity gradients also exert forces on illuminated ob- 
jects [8]. In this case, the dipole moment induced in the 
object responds to gradients in the field, yielding a force 
proportional to the gradient of the intensity. For a small 
sphere, this has the form [1, EH, 



F v (r) 



2 n 3 



k 2 a 



rrr 



VI. 



(12) 



Unlike the phase-gradient force, Fy can be directed up 
the optical axis, and so can be used to form stable three- 
dimensional optical traps, despite axial radiation pres- 
sure. This is the basis for optical tweezers and related 
single-beam optical traps [8] . 

Forces due to intensity gradients are manifestly con- 
servative [20[, arising as they do from the gradient of 
an analytic function. Phase-gradient forces, by contrast, 
generally are not conservative. This is evident because 



Vxg 



[in ri 



(VI) x (V$) + 0, 



(13) 



in general. This noteworthy property of scattering forces 
was pointed out by Ashkin [20[ , but appears not to have 
been emphasized subsequently in the optical trapping 
literature. Effects arising from non-conservative optical 
forces may not have drawn much attention because the 
experimentally accessible transverse forces can appear to 
be conservative in the symmetric beams that usually are 
considered. The phase-gradient force in a helical beam, 
by contrast, accelerates objects around the optical axis 
and so obviously does not conserve energy. 
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FIG. 1: Schematic representation of the optimized holo- 
graphic optical trapping system using shape phase holography 
to project extended optical traps. 



CREATING OPTICAL TRAPS WITH PHASE 
GRADIENTS 

We demonstrate and characterize phase-gradient forces 
using the optimized holographic optical trapping tech- 
nique 0, @, [2l[ and shape-phase holography [J 0] to 
create extended optical traps with specified profiles for 
both their intensity and phase. Our apparatus, shown 
schematically in Fig.[TJa), uses a phase-only spatial light 
modulator (Hamamatsu X7690-16) to imprint computer- 
generated holograms on the wavefronts of a slightly con- 
verging Q beam of light provided by a solid-state laser 
(Coherent Verdi 5W) operating at a vacuum wavelength 
of 532 nm. The modified beam is relayed to an objec- 
tive lens (Nikon Plan Apo, 100 x oil immersion, NA 1.4) 
that focuses it into the intended three-dimensional opti- 
cal trapping pattern. This objective lens is mounted in 
an inverted optical microscope (Nikon TE-2000U), and 
the trapping beam is directed into the objective's in- 
put pupil with a wavelength selective polarizing beam 
splitter. Light at other wavelengths passes through the 
beam splitter to form images on a CCD camera (NEC 
TI-324AII). 

The holograms used for this study are designed to bring 
the light to a focus along one-dimensional curves, C, em- 
bedded in the three-dimensional focal volume of the ob- 
jective lens. The hologram also encodes a designated 
intensity profile I(s) and phase profile (p(s) along the 
arclength s of C. This is accomplished by numerically 
back-projecting [22[ the intended field along C back to 
the plane of the SLM to obtain the ideal complex- valued 
hologram, \I/(r) = b(r) ex.p(ip(r)). The real- valued am- 
plitude b(r) is defined to be non-negative by absorbing 
sign changes into the definition of the phase, p(r). This 
facilitates encoding the complex hologram on a phase- 
only SLM 0, 0, H| • The shape-phase algorithm assigns 
p(r) to the phase pixel at r with a probability propor- 
tional to b(r). An alternate phase pattern imprinted on 
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FIG. 2: Linear dependence of a particle's mean velocity v 
along a uniformly bright line tweezer of length L on phase 
gradient q. 

imassigned pixels is used to divert excess light away from 
the intended trapping pattern. How the SLM plane is di- 
vided into assigned and unassigned domains can be tuned 
to minimize artifacts, improve diffraction efficiency, and 
optimize accuracy 0, 0, 0] • 

The images in Figs. 0Jb) and (c) show a focused line 
trap |2[ and ring trap, respectively, each designed to have 
uniform intensity and uniform phase gradients. These 
images were obtained by placing a mirror in the mi- 
croscope's focal plane and capturing an image of the 
reflected light. Because these extended traps come to 
a diffraction-limited focus, their axial intensity gradi- 
ents are steep enough to trap particles stably in three- 
dimensions against radiation pressure [8] . 



nal component of the radiation pressure. The measured 
angle of inclination and its uniformity along x confirm 
both the phase gradient's magnitude and also its unifor- 
mity. 

Line traps with tunable linear phase gradients were 
projected into aqueous dispersions of colloidal silica 
spheres 2a = 1.53 fim in diameter sealed into the 40 jam 
thick gap between a glass microscope slide and a #1 glass 
coverslip. The line traps were focused near the samples' 
midplane to minimize both interference with reflections 
from the glass-water interface and also hydrodynamic 
coupling to the walls. In that case, the spheres' viscous 
drag coefficient is reasonably approximated [24[ by the 
Stokes value 7 = 67rr]a, where 77 = 1 cP is the viscosity of 
water. Equations ([9]) and ([T0|) then suggest that a single 
sphere's speed, v, along the line should be proportional 
to q. 

To test this prediction, we measured the time required 
for a single sphere to travel the line's length, L, as the 
sign of q was flipped 20 times for each value of |g|. An es- 
timated 100 mW of light formed the trap, after account- 
ing for losses in the optical train. The observed off-line 
excursions of roughly 200 nm (root-mean-squared) sug- 
gest an axial and lateral trap stiffness comparable to that 
of a point-like optical tweezer powered by 1 mW. Under 
these conditions, the trapped sphere traveled the length 
of the line trap in roughly 2 s when subjected to the 
largest phase gradients. Results obtained by systemati- 
cally varying q are plotted in Fig. [21 and clearly show the 
anticipated linear dependence, except very near q = 0. 
The smallest values of |g| created too weak a force to 
overcome localized traps formed by uncorrected intensity 
variations. 



PHASE GRADIENT FORCES IN LINE TRAPS 
Linear phase gradients 

To demonstrate the phase gradient force predicted 
by Eq. ([9]), we used digital video microscopy to 
track micrometer-scale colloidal spheres moving along 
uniformly bright line traps such as that in Fig.[]Jb). The 
trap-forming hologram was adaptively optimized to min- 
imize intensity variations and to incorporate linear phase 
gradients, Vcp = qx, ranging from q = — 12 radians//im 
to q = 12 radians//im. The effect of imposing this trans- 
verse phase gradient can be seen in axial sections Q 
through the projected three-dimensional intensity dis- 
tributions, shown inset into Fig. O The diffraction- 
limited focal line remains in the xy plane despite the 
imposed phase gradient. The beam's direction of prop- 
agation, however, deviates from the z axis by the angle 
sin -1 (q/k). This tilt directs a component of the beam's 
wave vector along the line, thereby creating a longitudi- 



Phase-gradient traps 

More complicated phase gradients give rise to more in- 
teresting physical effects. The particles shown in Fig. [3] 
also are trapped along a uniformly bright line trap of 
length L = 10 /im. This line, however, has a parabolic 
phase profile along its length, ip(x) = ±(qx) 2 that is pre- 
dicted to create either force objects toward the center of 
the line, or out to the ends, depending on the choice of 
sign. The data in Fig. [3] also confirm this prediction. Ax- 
ial sections through the projected intensity pattern show 
that the phase-gradient barrier results from light diverg- 
ing along the line's length, while the well results from the 
projection of converging rays onto the trap's length. So 
long as the particles remain rigidly confined to the uni- 
formly bright focal line, the phase-gradient force approx- 
imates a conservative one-dimensional potential energy 
landscape. 
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optical vortex 



FIG. 3: Phase- gradient potential energy barrier (a) and well 
(b) created along a uniformly bright line trap. (1) Two 1.5 jam 
diameter silica spheres trapped on the line. (2) The uni- 
form in-plane intensity of the focused line. (3) Axial section 
through the measured intensity profile, show the divergence 
(a3) and convergence (b3) of the beams created by the im- 
posed phase profile. (4) The designed phase gradient, featur- 
ing the designed parabolic profile and off-line phase variations 
designed to minimize intensity variations. 



PHASE GRADIENT FORCES IN RING TRAPS 

Like holographic line traps, holographic ring traps, 
such as the example in Fig. [He), can be endowed with 
arbitrary phase profiles, including the uniform azimuthal 
phase gradient, ip(r) = £6, that defines a helical mode. 
A helical profile, by itself, causes a beam to focus into 
a ring of light, forming a torque-exerting optical trap 



)rque-exerti: 

known as an optical vortex 0, [H, Whereas the 

radius of an optical vortex is proportional to its helicity 
lolographic ring traps can be projected with any 
desired radius, independent of £. This facilitates system- 
atic studies of colloidal transport under varying phase 
gradients. Also unlike optical vortices, holographic ring 
traps come to a diffraction-limited focus, independent of 
i, and therefore can trap micrometer-scale objects stably 
in three dimensions. 

This property can be seen in the computed axial sec- 
tion through a holographic ring trap shown in Fig. Ufa). 
The section through the trap itself appears as two bright 
focal spots at the images mid-line. The ring's axial inten- 




-50 



1^ z t^ y 





-40 



-20 





t 



20 



40 



FIG. 4: Colloidal transport driven by azimuthal phase gra- 
dients in holographic ring traps, (a) Computed axial section 
through a holographic ring trap of radius R — 20 /am and he- 
licity £ = 30. (b) Volumetric representation of the measured 
three-dimensional intensity field in a holographic ring trap of 
radius R = 20 /im and £ = 10. Data points show the peak 
speed v of a single colloidal silica sphere circulating around 
the ring in (b) as a function of topological charge £. 



sity gradients are comparable to those of a conventional 
optical tweezer, and so suffice to trap objects stably in 
three dimensions. This ring trap features a helicity of 
£ = 30, and so is expected to exert a torque around 
the optical axis proportional to £. The helical wave- 
front topology also suppresses the beam's axial intensity 
through destructive interference. Consequently, none of 
the light projected by the objective lens appears along 
the beam's axis. Rather, it is diverted to a radius, Rg 
from the axis [H, [2tJ. Provided that the ring trap's ra- 
dius, i?, exceeds Rn, the converging helical beam focuses 
not only to the intended ring trap in the focal plane, 
but also to two conventional optical vortices above and 
below the ring trap. Sections through these optical vor- 
tices appear as bright regions above and below the mid- 
line in Fig. Ufa). Neither features the strong axial in- 
tensity gradients of the holographic ring. Consequently, 
neither can trap objects in three dimensions against ra- 
diation pressure. This structure also is evident in the 
three-dimensional reconstruction [3|] of an experimentally 
projected holographic ring trap in Fig. Efb). 

We measured the free-running speed of a single col- 
loidal silica sphere of radius a = 0.76 jam as it circulated 
around holographic ring traps of radius R = 2.6 /im pro- 
jected into the midplane of a 40 jam thick sample. The 
trapped particle was subjected to azimuthal phase gradi- 
ents ranging from £ = —50 to £ = 50, and its peak speed 
was measured [23| to within 10% for each value of the 
helicity. The results are plotted in Fig. 2J 
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Like optical vortices, holographic ring trap are sub- 
ject to £-fo\d and 2£-fold azimuthal intensity variations 
due to phase scaling errors [28]. These create local- 
ized traps that compete with the intended phase-gradient 
force [3, [2j|. Consequently, the particle remains mo- 
tionless for helicities below a threshold, t — i c . Similar 
pinning due to intensity gradients appears in Fig. O For 
t > £ C: however, the particle's peak speed increases lin- 
early with £ : consistent with the predictions of Eqs. ([9]) 
and ([TO]) . Intermittent circulation near t — i c gives rise 
to large velocity fluctuations characterized by giant en- 
hancement of the particle's effective diffusion coefficient 
[2^ . Disorder in the effective force landscape also gives 
rise to interesting collective dynamics for multiple par- 
ticles trapped on the ring, including transitions among 
periodic, chaotic and weakly chaotic steady states [301 ]. 
Phase-gradient forces in holographic ring traps therefore 
provide useful model systems for studying fundamental 
problems in nonequilibrium statistical mechanics. They 
also promise practical applications as the basis for mi- 
croscopic pumpsEH, mixers [32], and optomechanical 
micromachines 1331 . 



CONCLUSIONS 

The foregoing sections demonstrate that phase gradi- 
ents in beams of light give rise to transverse forces arising 
from the redirection of radiation pressure. The existence 
of orbital angular momentum in beams of light and its 
transfer to illuminated objects can be seen to be mani- 
festations of this general effect. 

Although phase-gradient forces may appear to be con- 
servative on restricted manifolds, they are not conser- 
vative in general. Non-conservative optical forces may 
give rise to interesting effects in illuminated particles' dy- 
namics, including departures from Boltzmann statistics 
for systems nominally in equilibrium. Such effects may 
substantially affect measurements of microscopic inter- 
actions and dynamics based on the statistics of optically 
trapped particles. They appear not to have been ad- 
dressed in the large and rapidly growing literature in this 
area. 

We have demonstrated that phase-gradient forces can 
be used to craft one-dimensional force profiles in uni- 
formly bright holographic line traps. These can be 
used to probe the interactions between multiple particles 
trapped on a line, including both intrinsic colloidal inter- 
actions and also light-induced interactions. This appli- 
cation for phase-gradient forces in holographic line traps 
will be discussed elsewhere. 

The phase-gradient force in holographic line traps is 
subject to the Abraham- Minkowski controversy. Ex- 
tended line traps created with shape-phase holography 
therefore may provide useful model systems for probing 
how a medium influences the momentum flux of light. 



Finally, it is worth mentioning that the polarization 
generally will vary with position in beams of light with 
spatially varying phase. Although optically isotropic 
materials are not influenced by polarization gradients, 
anisotropic materials can be. Polarization gradients 
therefore should provide additional independent avenues 
for exerting controlled forces on some microscopic sys- 
tems. 
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